Num         =   20;
PSI_0       =   1.2;
FracHN_w_0  =   1;
DiscFactor_0=   0.05;
FracTN_pC_0 =   1;
RES         =   SteadyState_Solver_Res(PP,SS,BellmanItSetup,[PSI_0;FracHN_w_0;DiscFactor_0;FracTN_pC_0]);
% Total Labor Supply
List_PSI    =   linspace(0.5,5,Num)';

List_RES    =   zeros(4,length(List_PSI));
parfor ii=1:Num
    List_RES(:,ii)  ...
                =   SteadyState_Solver_Res(PP,SS,BellmanItSetup,[List_PSI(ii);FracHN_w_0;DiscFactor_0;FracTN_pC_0]);
end

plot(List_PSI,List_RES(4,:));

% Total Labor Ratio
List_FracHN_w=  linspace(0.5,5,Num)';

List_RES    =   zeros(4,length(List_FracHN_w));
parfor ii=1:Num
    List_RES(:,ii)  ...
                =   SteadyState_Solver_Res(PP,SS,BellmanItSetup,[PSI_0;List_FracHN_w(ii);DiscFactor_0;FracTN_pC_0]);
end

plot(List_FracHN_w,List_RES(3,:));

% Saving
List_DiscFactor=linspace(0.05,0.2,Num)';

List_RES    =   zeros(4,length(List_DiscFactor));
parfor ii=1:Num
    List_RES(:,ii)  ...
                =   SteadyState_Solver_Res(PP,SS,BellmanItSetup,[PSI_0;FracHN_w_0;List_DiscFactor(ii);FracTN_pC_0]);
end

plot(List_DiscFactor,List_RES(1,:));

% Consumption
List_FracTN_pC=linspace(1,3,Num)';

List_RES    =   zeros(4,Num);
parfor ii=1:Num
    List_RES(:,ii)  ...
                =   SteadyState_Solver_Res(PP,SS,BellmanItSetup,[PSI_0;FracHN_w_0+1;DiscFactor_0;List_FracTN_pC(ii)]);
end
figure;
plot(List_FracTN_pC,List_RES(3:4,:));